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O ' Abstract 

O: 

We discuss a family of integrable systems on the sphere S with an 
additional integral of third order in momenta. This family contains the 
Coryachev-Chaplygin top, the Goryachev system, the system recently dis- 
covered by Dullin and Matveev and two new integrable systems. On the 
non-physical sphere with zero radius all these systems are isomorphic to 
£h ' each other. 



1 Introduction 

> . 

Let us consider particle moving on the sphere S 2 = {x G R 3 , \x\ = a}. As coordi- 
nates on the phase space T*S 2 we choose entries of the vector x = (xi,x 2 , x 3 ) and 
entries of the angular momentum vector J = p x x. The corresponding Poisson 
brackets reads as 



O 



{Ji j Jj } £-ijk'Jk i {Ji i 2>j / E-ijk%k j \ *&i j %j } , (1.1^ 



where is the totally skew-symmetric tensor. The Casimir functions of the 
j> \ brackets 

k> 1 3 3 

C3 ' i=l i=l 

are in the involution with any function on T*S 2 . So, for the Liouville integrability 
of the corresponding equations of motion it is enough to find only one additional 
integral of motion, which is functionally independent of the Hamiltonian H and 
the Casimir functions. 

If the corresponding Hamilton function H has a natural form, then according 
to Maupertuis's principle, integrable system on T*S 2 immediately gives a family 
of integrable geodesic on S 2 . If the additional integral of this integrable system is 
polynomial in momenta, integral of the geodesic are also polynomial of the same 
degree. 
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In this note we discuss a family of integrable systems on T*S 2 with a cubic 
additional integral of motion. Among such systems we distinguish the Goryachev- 
Chaplygin top PjJ |2| with the following integrals of motion 

H = J 2 + J\ + 4 J 3 2 + cxt , K = 2( J 2 + J 2 ) J 3 - CX3J1 , cel. (1.3) 

In j2j Chaplygin found the separated variables 

Qj = J 3 ± y/jfbJ$ + J$, J = 1,2, (1.4) 

in which dynamical equations are equal to 

(-iy ( qi - q 2 )q 3 = 2^P( qi ) 2 -a 2 c 2 q 2 , p(\) = A 3 - XH + K . (1.5) 

These equations are reduced to the Abel-Jacobi equations and, therefore, they 
are solved in quadratures 0. 

In variables qj (jl.4|) integrals of motion (jl.4j) read as 

H = q\ + q\q<i + q\ + cxi, K = qiq 2 (qi + 92) - cx 3 J\ . (1.6) 

In jH] Goryachev de facto substituted special generalizations of the variables qj 
(jl.4j) into expressions similar to (jl.6|) in order to construct new integrable system 
with a cubic integral of motion. In the next section we generalize this result. 

2 A family of integrable systems on the sphere 

Substituting canonical variables 

Qj = <*J3 ± ^Jf + Jl + f(x 3 )J 3 + g(x 3 ), {q u q 2 } = 0, (2.1) 
into the following ansatz for integrals of motion 

H = q\ + q x q 2 + q\ + m(x 3 )x 1 , 

(2.2) 

K = qiq 2 (qi + q 2 ) - n(x 3 ) J x - £(x 3 )xiJ 3 , 

one gets 

H = Jl + J 2 2 + (3a 2 + f{x 3 ))Jl + m{x 3 ) Xl + g(x 3 ) , (2.3) 

and 

K = -2aJ 3 (-a 2 J 2 3 + J 2 + J 2 2 + f(x 3 )J 2 + g{x 3 )} - n(x 3 ) J x - £(x 3 ) Xl J 3 . (2.4) 

Here a is an arbitrary numerical parameter, /, g, m, n and I are some functions 
of x 3 and of the single non-trivial Casimir a = y~x{ + x\ + X \ fll.2j) . 
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Theorem 1 On the phase space T*S 2 functions H \2. l J\) and K \2-4\) are in 

the involution with respect to the brackets if and only if function n(x 3 ) is 

solution of the following differential equation depending on a 2 

24a 2 - 9 15 Xsn ' - n " <K ° 2 + 3X3H " - n " ,(ya2 - X ^ (9 - nU " 



11 11' \ n' 2 

' (2-5) 

fhx*ri" - n""(a 2 - x 2 ) + 3n"\ 

+ n 



f5x 3 n'" -n""(a 2 - x 2 ) + 3n"\ 
V ) ' 



All another functions in \2.'J\2.4\) are parameterized by n(x 3 ) 



d n'(x 3 ) n(x 3 )n"(x 3 ) 

9{X3) = , x 2 , m{x 3 ) = , £{x 3 ) = — - — , 

n{x 3 y a n'(x 3 ) 

(2.6) 

2 3x 3 m(x 3 ) - 2(a 2 - x 2 )m'(x 3 ) x 3 £{x 3 ) - (a 2 - a%)£' (x 3 ) 



f(x 3 ) = 1 - 3a - a — r h 

n(x 3 ) n(x 3 ) 

Here d is arbitrary numerical parameter and z' = dz/dx 3 . 

The proof is straightforward. 

In this note we consider particular solutions of differential equation (j2.5j) only. 
Namely, substituting the following ansatz 

71(2:3) = c(x 3 + ef , c,e,/3eR, (2.7) 

in (|2.5|) one gets system of the algebraic equations on the three parameters a, (3 
and e whereas two other parameters c and d remain free. 

Theorem 2 Differential equation \2.5\) has five particular solutions in the form 
of P?7P only: 

1. ±a = (3 = 1, e = 0, ^(^3) = cx 3 , 

2. ±a = /? = -, e = 0, ^(^3) = ex]/ 3 , 

3. ±a = [3 = -, e = a, n(x 3 ) = c(x 3 + a) 1 / 6 , ^.8) 

4. ±a = /? = -, eGl, n(x 3 ) = c(x 3 + e) 1/2 , 

5. ±a = /? = -, e = a, n(x 3 ) = c{x 3 + a) 1//4 , . 
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The corresponding Hamilton functions are equal to 

d 



H 1 = J'( + J'2 + 4 J| + cx x 



1 ' 

X 3 



x 3 x 3 ' 



12 2(x 3 + a)J " (x 3 + a) s "> (x, + a) 1 ' 3 ' 



#4 = J I + J 2 + 1 + — " 77^ 7? ^3 + 7 ^777 + 

- 1 x 3 + e 4(x 3 + e) 2 / 3 (^s + e) 1 ^ X3 + e 



16 8(x 3 + a) / (x 3 + a) 3/4 (x 3 + a) 1/2 ' 

Transformation a — > —a /earfs to i/ie transformation of the free parameters 
(c,d) -> (-c, -d). 

Explicit expressions for additional cubic integrals of motion K\, . . . , K 5 may be 
obtained by using definition ()2.4|) and equations (j2.6l2.8j) . 

The Hamilton function Hi describes the Goryachev-Chaplygin top 2J. The 
second integrable system with Hamiltonian H 2 was found by Goryachev The 
Hamilton function H± and the corresponding cubic integral of motion K\ was 
studied by Dullin and Matveev [3]. The third and fifth integrable systems with 
Hamiltonians H 3 and H 5 are new. 

At present we do not know whether our systems in implicit or explicit forms 
(l2.5H2.9j) overlap with the families of integrable geodesic flows on S 2 considered 
by Selivanova [Hj and Kiyohara jSj. Recall that in j3HH] all the geodesic flows are 
defined in implicit form only (see also discussion in 0]). 



3 The Lax matrices 

In the fourth case (|2.8H2.9|) parameter e is free parameter and below we always 
put e = 0. 

Let us introduce 2x2 hermitian matrix 




where A is a spectral parameter and 

A(X) = {X-q l ){\-q2)=\*-2\aJ 3 +(a 2 -f{x 3 ))jl-Jl-J 2 2-g{x 3 ) ) 
B(X) = (x\ + ix 2 )m(x 3 )X + J 3 (xi + ix 2 )£(x 3 ) + (Ji + iJ2)n(x 3 ) , (3.1) 



D(X) = -n{x 3 ) 2 . 
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The trace of this matrix 

t(A) = A(X) + D(X) = X 2 - XH L + K L 
gives rise integrals of motion in the involution for the generalized Lagrange system 

H L = 2aJ 3 , K L = (a 2 - /(x 3 )) J 2 - J? - Jf - g(x 3 ) - n(x 3 ) 2 . 

The corresponding equations of motion may be rewritten in the form of the Lax 
triad 

jT(X)= [t(X),M(X)]+N(X), trJV(A) = 0. 

In contrast with the Lax pair equations at N(X) = 0, in the generic case deter- 
minant A(A) = detT(A) of the matrix T(A) (J3.1)) is dynamical function which do 
not commute with integrals of motion: 



1 

3 

1 
6 

1 
2 



A(A) 
A(A) 
A(A) 
A(A) 
A(A) 



a 2 / dn(x 3 ) Y 
P 2 V dx 3 J 



+ d, 



+ d, 



■0 (X + qi + q 2 ) 2 



-—X(X + q 1 + q 2 ) 



-- A(A + gi + q 2 ) 



dn 2 {x 3 ) 
dx 3 



+ d, 



(3.2) 



/ dn(x 3 ) \ ' 

V dx 3 J 



+ d, 



dn 2 (x 3 ) 
dx 3 



+ d. 



At ±a = (3=1 and n(x 3 ) = cx 3 matrix T(A) (j3.1j) was constructed in [Jj. In 
this case matrix T(A) defines representation of the Sklyanin algebra on the space 
T*S 2 associated with the symmetric Neumann system [7j. 

Theorem 3 If n(x 3 ) is one of the particular solutions 112. fy) of the differen- 
tial equations \2.5\i then T(A) \3. 1)) satisfies to the following deformation of the 
Sklyanin algebra 



{T(A), T(/*)} = [r(A - /*), T(A)r(/*) ] + Z(X, p) , 

1 2 

whereT(X) = T(X) ® / , = / ® T{^i), I is a unit matrix and 



(3.3) 



r(X-ii) 



2ia 
X — ji 



( 1 \ 

10 

10 

\ 1 / 



(3.4) 



Deformation Z(X,/a) is hermitian matrix 

( u(ji) 

U*(ji) 



Z(\,fi) 



\ 





-u*{\) w*(\,n) 




-u(X) \ 

w(X,fi) 



) 



(3.5) 



which depends of the entries ofT(X) h3. 1)) only: 
• at (3 = 1 we have u = w = 0; 
1 1 



at (3 



at p 



3' 6 
u(ji) = 

1 1 



we have 



-Aia 



\fWF 7 d \/ A(/x) - dB(\) - a/ A(A) - dB(fi) 



-Aid 



A(A) - A(ji) 
A — ji 



A — [i 



2' 4 



we have 



-2ia 



A(/i) - dfiB(X) -XB(fx) 
/iD(fj,) X — /i 

A(A) - A(/i) 



-2za 



A — n 



Here A(A) = A(A)A(A) - B(X)B*(X) is determinant of the matrix T(A) / TO) . 

The proof is straightforward. 

One of the main properties of the Sklyanin algebra is that for any numerical 
matrices K, and for some special dynamical matrices K, coefficients of the trace of 
matrix J£?(A) = fCT(X) give rise the commutative subalgebra 

{tr/CT(A) ,tr/CT(/i)} = 0, 

(see jH] and references within). All the generators of this subalgebra are linear 
polynomials on coefficients of entries Tij(X), which are interpreted as integrals of 
motion for integrable system associated with matrices T(A) and K, jH]. 

Deformation of the Sklyanin algebra ()3.3t 13. 5jl has the similar property. 



Theorem 4 // dynamical matrix K, has the form 



X + 2aJ 3 6i 
ci 



bi, Ci G C 



then coefficients of the polynomial 

P(X) =tr/CT(A) = X 3 -XH + K 
are in the involution on T*S 2 . 



(3.6) 
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If bi — Ci — 1/2 then the first coefficient H in ^(A) (|3.6|) coincides with one of 
the Hamiltonians Hi, . . . , H5 ()2.9|) listed in the Theorem 2, whereas the second 
coefficients K is the corresponding cubic integral Ki, . . . ,K 5 ()2.4|) . If bi and c\ 
is arbitrary one gets the same Hamiltonians up to the suitable rescaling of x and 
rotations 

x^bUx, J^UJ, (3.7) 

where b is numerical parameter and U is orthogonal constant matrix. 

The equations of motion associated with the Hamilton function H (j3.6|) may 
be rewritten as a Lax triad for the matrix S£[X) = /CT(A) 

= + tr^(A) = . 

Here matrices ^ and jY are restored from the deformed algebra (J3.3|) and defi- 
nition of Hamiltonian (|3.6|) in just the same way as for the usual Sklyanin algebra 



4 Isomorphism of the systems at a = 

For all the considered systems (|2.9|) at a = additional term Z(\,fi) in (|3.3|) is 
equal to zero according to ()3.2I3.5|) 

a = =^ A(X)=d => Z(A,aO = 0. 

In this case matrices T(A) associated with five integrable systems ()2.9|) define 
five representations of the Sklyanin algebra on the space T*S 2 . Of course, these 
representations are related to each other by canonical transformations. 

Theorem 5 At a = 0, i.e. on the non-physical sphere S 2 with zero radius, 
integrable systems listed in the Theorem 2 are isomorphic to each other. 

To prove this Theorem we introduce variables 

Pi = ^-hii%)= (4.8) 

= — \n^qj(xi +ix 2 )m(x 3 ) + J 3 (xi + ix 2 )£{x 3 ) + (Ji + iJ2)n(x 3 )j . 
At a = variables pi >2 and q± t 2 are canonical Darboux variables according to ()3.3|) 

{Pi, Qj} = $ij , {Pi, Pj} = {<li,qj} = Q, h 3 = 1, 2. 

In order to construct canonical transformations which relate integrable systems 
with Hamiltonians Hi, . . . , H 5 ()2.9|) we have to identify variables pi j2 , qi,2 ()2.1I4.8|) 
associated with the different functions n(x 3 ) (|2.8j) . 
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We could not lift these symplectic transformations to the Poisson maps. So, 
we can not assert that integrable systems ()2.9|) are isomorphic on the generic 
symplectic leaves (II. 2j) . 

The result of the Theorem 5 may be interpreted in the following way. At a = 
on the special symplectic leaf of the Lie algebra e(3) there exists a germ of single 
integrable system with Hamiltonian H (J1.3|) . Using canonical symplectic trans- 
formations one can gets infinitely many different forms of this integrable system. 
However, according to the Theorem 2, these different forms of the germ admit 
only a denumerable set of the continuation on the generic symplectic leaves with 
conservation of the integrability property. The similar observation for another 
family of integrable systems on the sphere is discussed in 

5 Summary 

Using the separation of variables for the Goryachev-Chaplygin top we constructed 
a family of integrable systems on the sphere with a cubic additional integral 
of motion. On the non-physical sphere S 2 with zero radius these systems are 
isomorphic to each other. 

On this non-physical sphere a = the separated variables for all five sys- 
tems coincide with the separated variables for the Goryachev-Chaplygin top up 
to symplectic transformations. It allows us to integrate equations of motion in 
quadratures. On the usual sphere at a ^ the separated variables are unknown. 
We suppose that these variables may be constructed using the proposed defor- 
mation of the Sklyanin algebra. 

The research was partially supported by RFBR grant 02-01-00888. 
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